Ergodicity properties of quantum expectation values in entangled states 
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Using a model Hamiltonian for a single-mode electromagnetic field interacting with a nonlinear 
medium, we show that quantum expectation values of subsystem observables can exhibit remarkably 
diverse ergodic properties even when the dynamics of the total system is regular. The time series 
of the mean photon number is studied over a range of values of the ratio of the strength 7 of the 
nonlinearity to that of the inter-mode coupling g. We obtain the power spectrum, estimate the 
embedding dimension of the reconstructed phase space and the maximal Liapunov exponent A max , 
and find the recurrence-time distribution of the coarse-grained dynamics. The dynamical behavior 
ranges from quasiperiodicity (for j/g <C 1) to chaos as characterized by A max > (for j/g > 1), 
and is interpreted. 
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Classical nonlinear dynamical systems generically ex- 
hibit a long-time behavior that is ergodic — either on the 
energy surface (in Hamiltonian systems) or on some at- 
tractor (in dissipative systems). A hierarchy of random- 
ness can be identified, ranging from ergodicity to global 
exponential sensitivity to initial conditions. Quantum 
systems, on the other hand, are governed by linear equa- 
tions (the Schrodinger equation for the state vector or 
the Liouville equation for the density operator). The 
randomness in a quantum system, and the precise man- 
ner in which its information content changes with time, 
are not easily determined. 

At least two different approaches have been employed 
in identifying signatures of the ergodicity properties of 
the quantum counterpart of a generic classical system 
[H . One of these relies on the observation that such sig- 
natures are manifested in the energy-level statistics of the 
corresponding quantum system. If the system is clas- 
sically integrable, the quantum levels cluster together, 
and could even cross when a specific parameter in the 
Hamiltonian is varied Q. A classically chaotic system, 
on the other hand, has its corresponding quantum levels 
so correlated as to resist such crossings [3(. Another ap- 
proach is based on the dynamics of the overlap between 
two quantum states of the same physical system which 
originate from the same initial state, but with slightly 
different values of one of the control parameters [J] . The 
time-dependent overlap is close to unity for all t if the 
normalized initial state is located in a regular region of 
the classical phase space. In contrast, if the initial state 
is in a chaotic region of this space, the overlap falls off 
exponentially in time. 

These lines of investigation concern quantum signa- 
tures of classical dynamics. The inverse problem is of 
importance, and has also received attention: namely, the 
identification of signatures of non-classical effects such as 



wave packet revivals in the temporal behavior of quantum 
expectation values (which, in turn, could be regarded 
as effective dynamical variables in an appropriate 'phase 
space'). (A (near-)revival of an initial state | -0 (0) ) at 
t = T rcv implies that | (^>(0)|^(T rev )) | 2 ~ I.) The dy- 
namics of a single mode of the radiation field governed 
by a nonlinear quantum Hamiltonian H enables us to un- 
derstand the connections between the behavior of quan- 
tum expectation values and various non-classical effects 
displayed in wave packet dynamics 5]. 

In order to explore more thoroughly the range of dy- 
namical behavior of expectation values of observables, in 
particular in the presence of entanglement, we need a sys- 
tem in which revival phenomena can either occur or be 
suppressed, depending on the values of the parameters 
in H. An uncomplicated but nontrivial H for our pur- 
poses is the one that describes the interaction of a single- 
mode field of frequency u with the atoms of the nonlin- 
ear medium through which it propagates. The medium is 
modeled [6[ by an anharmonic oscillator with frequency 
ujq and nonlinearity parameter 7. The Hamiltonian of 
the total system is given by 

H = hwa)a + huj Q tfb + fir/ 6 t2 b 2 + hg (a f 6 + 6 f a). (1) 

(a,a + ) are the field annihilation and creation operators, 
(fo, 6+) are the corresponding atomic oscillator operators, 
and g quantifies the coupling between the field and atom 
modes. Importantly, although the photon number oper- 
ator N = a^a does not commute with H for any g 7^ 0, 
the total number operator TV* ' = (et + a + b^b) does so 
for all values of the parameters in H. While this implies 
that H can be cast in block-diagonal form in a direct- 
product basis of field and atom Fock states, the model is 
not trivial. A simple way to see this is to re-write H in 
terms of spin operators Ji (i — 1,2,3) constructed from 
the two pairs of boson operators (a, a^) and (6, b^). It is 
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then easily seen that the system of Heisenberg equations 
of motion for the spin operators does not close. 

The quantum mechanical expectation value (N(t)) (or 
the mean energy of the field mode) serves as a very conve- 
nient variable to probe the dynamics of this subsystem. It 
varies with time because of the coupling between the two 
modes, and it deviates from periodicity because of the 
nonlinearity in H . Its temporal behavior is remarkably 
diverse, ranging from quasiperiodic to chaotic, depending 
strongly on the initial state and the parameter regime — 
in particular, on the degree of coherence of the initial 
state, and on the ratio 7/g of the strengths of the nonlin- 
earity and the field-atom coupling. We have studied the 
dynamics for states that are initially non-entangled direct 
products of the field and atomic oscillator states: specifi- 
cally, states with the field in a coherent state \a) (CS) or 
an m-photon-added coherent state \a,m) (PACS), while 
the atomic oscillator is in its ground state |0). Recall 
that the CS |a) (a S C) satisfies a \a) = a\a), and is 
a minimum uncertainty state. The normalized PACS is 
defined Q as \a,m) = [a)) m \a) /[to! L^-v)] 1 / 2 where 
to is a positive integer, v = |a| 2 , and L m is the La- 
guerre polynomial of order m. A PACS possesses the 
useful properties of quantifiable and tunable degrees of 
departure from perfect coherence and Poissonian photon 
statistics. For brevity, we write \a) (g> |0) = \a; 0) and 
\a,m) (g) |0) = I (a, to); 0) for the initial states consid- 
ered. (N(0)) = v and [(to + 1) L m+1 {-v) / L m {-v)\ - 1, 
respectively, in the two cases. 

We have carried out a detailed analysis of the time se- 
ries (using a time step 8t ranging from 10 -2 for small 
7/5 to 10 _1 for large 7/3) generated by the values of the 
mean photon number (N(t)) computed over long inter- 
vals of time (10 6 time steps), including phase space recon- 
struction, estimation of the minimum embedding dimen- 
sion demb, calculation of the power spectrum 
and recurrence-time statistics (using time series of 10 7 
steps when necessary). We use a robust algorithm [ll| 
developed for the estimation of the maximal Liapunov 
exponent A max from data sets represented by time series 
(see also [iH). The phase-space reconstruction procedure 
(including the extraction of d em b) has been carried out 
carefully, and it has been checked that any further in- 
crease in the dimensionality of the reconstructed phase 
space does not alter the inferences made regarding the 
exponential instability, if any, of the system. 

For small values (<§; 1) of 7/5, near-revivals and frac- 
tional revivals of the initial state occur, that are mani- 
fested in the entropy of entanglement of the system [3] • 
Correspondingly, we find that the dynamics of (N(t)) 
ranges from periodicity to ergodicity, but is not chaotic, 
essentially independent of the nature of | -0 (0)) . As a 
case representative of weak nonlinearity, we have chosen 
the parameter values 7 = 1, g = 100. Figure [T] (a), a 
log-linear plot of the power spectrum S(f) (the Fourier 
transform of the autocorrelation computed from the time 



series of (N(t))) as a function of the frequency / when 
\ip(0)) = \a; 0), is indicative of quasiperiodic behavior. 
With increasing lack of coherence of the initial field state, 
the number of frequencies seen in £>(/) increases. This 
is already evident from Fig. [T] (b), which corresponds to 
M0)> = |(«>5);0>. 
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FIG. 1: Power spectrum of the mean photon number vs. the 
frequency (in units of g) for the initial states (a) \a; 0) and 
(b) I (a, 5) ; 0) with j/g = 10 -2 and v = 1. 

In contrast to the case of weak nonlinearity, the na- 
ture of the subsystem dynamics changes drastically when 
l/9 ~ !■ A s representative values for this nonlinearity- 
dominated regime, we have set 7 = 5,5 = 1. We first 
examine the case corresponding to v = 1. For an ini- 
tial field CS, both the time series and S(f) confirm that 
the subsystem dynamics is not chaotic. In contrast to 
this, an initial PACS leads to a chaotic form for S(f), for 
sufficiently large values of to. This is supported by an es- 
timation of A max from the time series. The initial set of 
separations between the j th pair of nearest neighbors in 
the reconstructed phase space evolves to the set {dj(k)} 
after k time steps. Then A max is the slope of the plot 
of (In dj(k)) (the average is over all values of j) against 
t in the linear region lying in between the initial tran- 
sient and final saturation regions. Figure [2] (a) depicts 
(lndj(fc)) vs. t for IV'(O)) = |(a, 5) ; 0), and the estimate 
of A max , whose positivity indicates a chaotic variation of 
the mean energy of the field mode. For a given value of 
"f/g, an increase in v leads to chaotic behavior even for an 
initial coherent field state. This is demonstrated in Fig. 
0(b), which shows (Indjik)) vs. t for |^>(0)) = \a; 0) 
with v = 10. We further find that A max increases sys- 
tematically with to for an initial PACS. 

Table 1 summarizes these conclusions. In order to rule 
out round-off or truncation errors as the source of the 
computed chaotic behavior, we have verified in each case 
that the conlcusions are not altered if {b'b} is chosen as 
the signal for which the time-series data is computed, and 
that {N} + (b^b) does remain constant, as required. The 
entropy of entanglement (not presented here) provides 
independent corroboration of the dynamical behavior as 
deduced from S(f) and A max . 

Our conclusions are reinforced by a detailed analysis 
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I V(0)> = |(a,l); 0), v = 1. According to Table 1, this 
case is non-chaotic. Figure|3](a) shows the invariant den- 
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FIG. 2: (In dj(k)) vs. t for the initial states (a) |(q,5) ; 0) 
with j/g = 5 and v = 1 and (b) \a ; 0) with j/g = 5 and v — 
10. The solid line corresponds to an embedding dimension 
dcmb = 5, and the dotted lines to values of d om b from 6 to 10. 
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Table 1: Qualitative dynamical behavior of the mean photon 
number of a single-mode electromagnetic field interacting 
with a nonlinear medium. "Regular" =>• A max = 0. 



of another important characterizer of dynamical behav- 
ior: recurrence statistics of the coarse-grained dynam- 
ics of (N(t)) as represented by its time series. For the 
range of parameter values we use, the scatter in (N(t)) 
is typically > 1. We use a cell size ~ 10 -2 and very 
long time series. This enables us to numerically con- 
struct the invariant density p (and hence the stationary 
measure /i for any cell C), as well as the distribution 
F(t) of the time t of first recurrence or return to C. 
The mean recurrence time (r) can then be calculated, 
and compared with the result (r) = that follows 
from the Poincare recurrence theorem. 14j. As the lat- 
ter is derived from the requirement of ergodicity alone, 
an agreement between the two values confirms that the 
dynamics is indeed ergodic in all the cases studied. We 
present here just two representative cases, both of which 
are also included in Table 1, for ready reference. The 
first corresponds to weak nonlinearity (j/g = 1CP 2 ) and 
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FIG. 3: (a) Invariant density and (b) first-recurrence-time 
distribution for the cell C from the time series of {N(t)), for 
weak nonlinearity. F(t) is characteristic of quasiperiodic dy- 
namics. 



sity, while (b) shows the actual recurrence time distribu- 
tion. The discrete nature of the latter is a clear indica- 
tion that the dynamics is actually quasiperiodic [l5j |. In 
marked contrast, consider a case of strong nonlinearity, 
7/5 = 5 and |-0(O)) = \(a, 1) ; 0) , v = 10. According to 
Table 1, this case is chaotic, with A max = 0.80. Figure 



C=[l 1.39,11.41] 





FIG. 4: (a) Invariant density and (b) first-recurrence-time 
distribution for the cell C from the time series of {N(t)), for 
strong nonlinearity. F(r) is characteristic of chaotic dynam- 
ics. 



|4] (a) shows that the invariant density is in fact well- 
approximated by a Gaussian in this case. More impor- 
tantly, F(t) is very well fitted by the exponential distri- 
bution fi e _MT . This is precisely the distribution expected 
in a hyperbolic dynamical system, for a sufficiently small 
cell size [lfl]. Moreover, in such a system successive re- 
currences to a cell must be uncorrelated and Poisson- 
distributed. We have further confirmed this feature by 
examining the distribution of two successive recurrences, 
using even longer time series (10 7 steps). The distribu- 
tion is again well fitted by the next term in the Poisson 
distribution, /i 2 re _A1T . 

Now consider the completely classical counterpart of 
the Hamiltonian in Eq. ([T]). Let the linear harmonic 
oscillator associated with (a, a') have a mass m, position 
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x and momentum p x , and let that associated with (b, tf) 
have a mass M, position y and momentum p y . Putting 
in all the constant factors (including h) in the definitions 
of the raising and lowering operators in Eq. (fTJ) , we get 
H{x, p x , y, p y ). When h — > 0, the only consistent way 
to obtain a non-trivial, finite expression for the classical 
Hamiltonian H c \ is to let 7 — > simultaneously, such that 
the ratio j/H — > A = a finite number. Then, with H± = 
p 2 x /(2m) + muj 2 x 2 /2 and H 2 = p 2 y /{2M) + Mu 2 ,y 2 /2, we 
find 

flci = H 1 + H 2 + (X/uj 2 )Hi + 

+ (g/^/uJuJo) rriM ujluq xy + p x Py/VmM^J . (2) 

The counterpart of N tot is = Hi/u + H 2 /ujo , 

which Poisson-commutes with H c \ . Hence the 2-freedom 
classical system is Liouville- Arnold integrable. Further, 
although H c \ has cross terms that could change sign, 
7V'[ 0t = constant is a hyperellipsoid in the 4-dimensional 
phase space, so that the motion is bounded for any set of 
initial conditions. All four Liapunov exponents vanish, 
and the classical motion is always regular, and restricted 
to a 2-torus, for each set of initial conditions. This behav- 
ior is indeed very different from the much more diverse 
one found for the quantum expectation value (N(t)}. 

What, then, is the interpretation of a positive value for 
A max , and the hyperbolicity implied by the recurrence- 
time statistics, as deduced from the time series for the 
subsystem variable (N) (equivalently, (b^b)) in the quan- 
tum mechanical case? To start with, we note that a 
comparison of the quantum and classical cases is not al- 
ways straightforward [rTj . and the case at hand is one 
such instance. Letting 7^0 would remove the nonlin- 
ear term b' 2 b 2 in the quantum H . In that case the dy- 
namics of (N(t)) reduces to a trivial periodic exchange 
of energy between the field and atom modes. The ori- 
gin of this dichotomy can be traced back to the inad- 
equacy of the naive Ehrenfest theorem, which does not 
generally take into account the non-commutativity be- 
tween x and p, in retrieving the classical regime of the 
quantum system. An outcome of this feature is that the 
Liapunov exponents that characterize the dynamical be- 
havior of classical and quantum expectation values of the 
same observable can indeed be very different from each 
other 18]. The following observation [l9[ is pertinent 
in this regard. In isolated quantum systems with a dis- 
crete energy spectrum, using unitarity and the Schwarz 
inequality, it can be established that the Liapunov ex- 
ponents would vanish, if computed from time-series data 
collected over a sufficiently long time (which could be 
much longer than the characteristic time scales in the 
problem), indicative of non-chaotic behavior. However, 
once measurement upon the system is included through 
appropriate interaction with an external system, the cor- 
responding Liapunov exponents need not vanish. In the 
present context, the interaction between the modes is ef- 



fectively tantamount to continual measurement upon ei- 
ther subsystem. Hence the dynamics of a subsystem, as 
deduced from time-series data of the corresponding vari- 
ables, may show chaotic behavior, even if the system as 
a whole does not. The exponential instability associated 
with a positive Liapunov exponent is indicative of the 
manner in which an initial wave packet spreads, and the 
entanglement of the system increases, under time evolu- 
tion. These aspects are worth bearing in mind, inasmuch 
as systems are ultimately quantum mechanical, and mea- 
sured data are generically the time series of observables. 
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